Hydrodynamic instabilities provide a generic route to 
spontaneous biomimetic oscillations in active filaments 
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Microscopic living organisms locomote in viscous fluids through spontaneous beating of filamen- 
tous structures anchored at one end, such as flagella and cilia. Prokaryotic flagella rotate rigidly 
like a corkscrew, while eukaryotic flagella are flexible and oscillate in a plane. We observe similar 
biomimetic beating behavior in silico by clamping an active filament at one end and solving for hy- 
drodynamics interactions. Highly active filaments become unstable to transverse perturbations and 
exhibit autonomous oscillations. This transition into a limit cycle occurs via a supercritical Hopf 
bifurcation. The time period and amplitude of beating increase with increasing filament length, but 
collapse to a master curve on appropriate scaling. An analytical calculation of the spectrum of the 
filament model in the free-draining approximation fails to capture oscillatory behavior, emphasizing 
the crucial role played by hydrodynamic interactions. 



Nonequilibrium processes by which biochemical energy 
is converted to mechanical stress are essential for life. 
These processes enable the autonomous motility of liv- 
ing organisms which, otherwise, would be impossible in 
equilibrium. At the microscale the engines of motility 
are primarily cilia and flagella, which are bundles of in- 
extensible elastic filaments that oscillate when forced by 
the energy transduction of molecular motors. Although 
the mechanism behind this behavior has been studied in 
varying degrees of detail over the years [1-12], the contri- 
bution of hydrodynamic interactions (HI) has generally 
been overlooked. The role played by HI in biomimetic 
oscillations in synthetic active filaments has been stud- 
ied in recent experiments [13, 14]. Motivated by these 
biological and biomimetic examples we seek to theoret- 
ically understand the contribution of HI to autonomous 
motion in active filaments. Here, we show that the 
competition between the destabilising effect of nonequi- 
librium hydrodynamic flow and the stabilising effect of 
bending and stretching elasticity provides a generic route 
to spontaneous biomimetic beating in a minimal model 
of active filaments [15]. Our results show that HI due 
to force-free and torque-free hydrodynamic flows gener- 
ated by active energy transduction is sufficient to induce 
and sustain spontaneous oscillations in an inextensible 
filament. This minimal physical requirement then im- 
plies that biomimetic beating, previously observed only 
in bundles of active filaments (BAMs) [13, 14], can be 
replicated in a single filament, for instance, a chain of 
autocatalytic nanorods [16-19] or other synthetic active 
particles. Our results will help answer scientific ques- 
tions like what minimal physical ingredients are required 
to generate synchronised metachronal waves, help guide 
the rational technological design of nanoscale machines 
like carpets of beating filaments that enhance fluid mix- 
ing, and bring us closer to uncovering the nonequilibrium 
mechanisms that drive biological motility. 

Prokaryotic bacteria [20] as well as eukaryotic sperm 
cells [21, 22] employ rhythmic flagellar beating for loco- 



motion in viscous fluids. Bacterial flagella rotate rigidly 
in corkscrew fashion [23, 24], whereas spermatic flag- 
ella behave more like flexible oars [25] with their beating 
mostly confined to a plane [26, 27]. Oscillatory motility 
in clamped flagella often arise spontaneously and, with an 
unlimited supply of energy, can persist indefinitely with- 
out any external or internal regulatory pacemaker mech- 
anism [22, 28]. Autonomous motility as well as sponta- 
neous beating due to hydrodynamic instabilities has been 
recently reproduced in vitro [13, 14], where a biomimetic 
active motor-microtubule assemblage has been shown to 
exhibit remarkable cilialike beating motion with hydro- 
dynamic interactions playing a crucial role in synchro- 
nized oscillations [13]. 

In this Letter, we study a minimal active filament 
model [15] which, once clamped at one end, exhibits a 
variety of spontaneous beating phenomena in a three di- 
mensional fluid. Our model filament consists of active 
monomeric beads in a viscous fluid connected through 
spring potentials that enforce inextensibililty, semiflexi- 
bility and self- avoidance. Activity is modelled through 
force- free and torque- free singularities [29-34], since the 
net force and torque exerted on the fluid by any inter- 
nally active system must vanish. The most dominant 
of these, the stresslet, is retained in our study. Hydro- 
dynamic interactions are propagated by the appropriate 
Green's functions. The equation of motion for the ac- 
tive filament [15] incorporating the effects of nonlinear 
elastic deformations, active processes and hydrodynamic 
interactions is 

N 

Tn=Yl r ™) ' f ™ + D ( rn r ™) ' a ™] (1) 

m—1 

where r n is the location of the n-th bead, f n is the total 
elastic force on the n-th bead, and cr n = cr (t n t n — 1/3) is 
stresslet tensor directed along the the local unit tangent 
t n . Here <r sets the scale of the activity. The monopolar 
Oseen tensor O and the dipolar stresslet tensor D re- 
spectively propagate the elastic and active contributions 
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FIG. 1. (Color) Biomimetic oscillations of the clamped filament plotted at different times over an oscillation period T. In (a) 
we see rigid aplanar corkscrew rotation for A — 100 while in (b) we see flexible planar beating for A — 200. The color of the 
beads as well as the trace of the tip correspond to individual instantaneous monomer speeds. 



to the flow. We impose clamped boundary conditions at 
one end and solve the equation of motion through direct 
summation of the hydrodynamic Green's functions. For 
a filament of length L and bending modulus n the dy- 
namics is characterised by the dimensionless measure of 
activity A = Lao/ft, called the activity number [15]. 

Here we show, using linear stability analysis and nu- 
merical simulations, that the active filament has a hy- 
drodynamic instability with increasing activity number. 
The containment of this linear hydrodynamic instability 
by the non-linear elasticity leads to spontaneously oscil- 
lating states. We briefly recall the mechanism of this 
linear instability [15]. Extensile activity produces flows 
that point outward in the longitudinal direction and in- 
ward in the normal direction. A spontaneous transverse 
perturbation breaks flow symmetry about the filament 
midpoint, resulting in a net flow along the outward fila- 
ment normal. At low A, this is countered by the elastic 
restoring force. Beyond a critical activity number A&, 
linear elasticity is unable to balance the outward flow 
leading to a linear instability. The instability produces 
large filament deformations which are ultimately con- 
tained by the non-linear elasticity. Constrained by the 
clamp, the filament acquires a rigid corkscrew conforma- 
tion and rotates about the initial axis. This rotational 
corkscrew motion is reminiscent of prokaryotic flagellar 
beating [23, 24]. We show this motion in Fig. 1(a) over 
one time period of oscillation together with the projec- 
tion of the filament on the plane of the clamp. A section 
of the three-dimensional flow in a plane perpendicular to 
the plane of the clamp is shown in Figs. 2(a) and 2(b). 
The net flow points in the direction opposite to the fil- 
ament curvature and the entire flow pattern co-rotates 



with the filament. 

With increasing A the filament acquires greater flex- 
ibility. At a second critical value of the activity A C 2, 
the filament motion collapses to a two-dimensional plane 
perpendicular to the plane of the clamp in which it beats 
spontaneously. This symmetric beating motion is rem- 
iniscent of eukaryotic flagellar beating [22, 26-28]. We 
show this motion in Fig. 1(b) over one time period of 
oscillation together with the projection of the filament 
on the plane of the clamp. The projection is now a line, 
showing that motion is confined to a plane. A section 
of the three-dimensional flow in the plane of beating is 
shown in Figs. 2(c) and 2(d). Two distinct types of fil- 
ament conformations of opposite symmetry are now ob- 
served. In the even conformation shown in Fig. 2(c), the 
flow points in the direction opposite to the curvature as 
in the corkscrew state. However, in the odd conforma- 
tion shown in Figs. 2(d), the flow has a center of vorticity 
at the point of inflection of the filament. This center of 
vorticity moves up the filament and is shed at the tip at 
the end of every half cycle. 

In Fig. 3(a) we show the variation of the time period 
of oscillation T for both the corkscrew rotating and pla- 
nar beating states for a filament of fixed length L. While 
the period T depends jointly on n and <Jo, a collapse is 
obtained when the period is scaled by the active relax- 
ation rate T a = ao/r]L 3 , where n is the viscosity of the 
fluid. The oscillations become more rapid with increasing 
activity. In Fig. 3(b) we show the variation of the real 
parts of the two largest eigenvalue pairs obtained from 
a linear stability analysis. Here too, the data collapses 
when scaled by the active relaxation time T" 1 . All eigen- 
values are real and negative as A is increased from zero. 
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With increasing A, the two largest eigenvalue pairs ap- 
proach, converge, and become complex conjugate pairs. 
This corresponds to a transition from a stable node to a 
stable focus where the response changes from being over- 
damped to underdamped. With further increase of A the 
complex eigenvalues approach the imaginary axis mono- 
tonically, crossing them at a critical value A c i- Through 
this supercritical Hopf bifurcation, the stable focus flows 
into the limit cycle corresponding to the corkscrew rota- 
tion. The variation of the real parts of the first eigenvalue 
pair reveals a striking feature of this filament. Since ex- 
tensile activity produces flows that enhance the pertur- 
bation, the real parts for the most part increase with 
A, according to expectation [11]. However, for a small 
but significant range of values of A, we find that this 
quantity decreases with activity. Thus, in this window, 
extensile activity suppresses perturbations, which thus 
relax faster. With further increase of A beyond A C 2 the 
corkscrew limit cycle transforms into the limit cycle cor- 



responding to the planar beating state. Our simulations 
provides evidence to support mult ist ability and hystere- 
sis in this transition, the precise nature of which will be 
quantified from a Floquet stability analysis in a future 
study. 

The phenomena described above remain qualitatively 
the same as the length of the filament is changed. How- 
ever, both the time period and amplitude of the oscilla- 
tions increase with the filament length, as shown in insets 
of Figs. 4(a) and 4(b). Remarkably, both these quan- 
tities appear to obey a scaling relation L~ a f(A/L@), 
as shown in the main panels. We estimate the expo- 
nents to be a = —f3 = 1.3 for the time period and 
a — —1.5,13 = —1.3 for the amplitude. 

To better understand the role of hydrodynamic inter- 
actions, we have performed a linear stability analysis on 
the free-draining continuum limit of this model. Our re- 
sults, obtained analytically, show that the spectrum of 
the clamped filament system in this limit is purely real, 



4 




100 150 

A 

(a) 



juiescent 



Stable node , Stable focus 



Beating 
Limit cycle 



20 



40 

A 



(b) 



A 


a = 0.2 





a = 0.4" 


o 


k = 0.5 


a 


K= 1.0 



60 



80 



FIG. 3. (Color online) Variation of (a) the scaled timeperiod of filament beating and (b) the real parts of its primary eigenvalues 
with the activity number A, plotted for various values of k and ao. In (a) we show the appearence of spontaneous oscillations in 
the filament at A ~ 50 corresponding to rigid corkscrew rotation, followed by a transition at A ~ 180 to flexible planar beating. 
In (b) we demarcate the various dynamical regimes of the filament obtained from linear stability analysis. The collapse of the 
different plots onto a single master curve for both (a) and (b) shows that A completely determines filament dynamics at fixed 
L. 




FIG. 4. (Color online) Variation of the rescaled (a) timeperiod and (b) amplitude of filament beating with the length L and 
activity number A. The unsealed results in the insets show that the time period as well as the amplitude of beating increases 
with L, while the rescaled main panels show that these variations are well captured by a scaling form L~ a f(A/ L^). Our 
estimates are a — —f3 — 1.3 for the time period and a — — 1.5,/? = —1.3 for the amplitude. 



with the transition to instability occurring via a saddle- 
node bifurcation, as seen in Fig. (5). This rules out any 
oscillatory modes of the filament. Since the free-draining 
approximation only takes local hydrodynamics into ac- 
count, this analysis suggests that hydrodynamic interac- 
tions are required to produce an oscillatory instability in 
our model active filament. 

It is thus quite extraordinary that this minimal 



nonequilibrium active filament model shows spontaneous 
emergent biomimetic behavior which is reminiscent of the 
rhythmic oscillations of various prokaryotic and eukary- 
otic flagella [20-27]. The generic nature of our model 
indicates that chains of elements which convert chemical 
energy to mechanical stresses in a fluid environment, for 
instance synthetic catalytic nanorods [16-19], should ex- 
hibit the behaviour described here. Our results lead us to 
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conclude that hydro dynamic instabilities due to internal 
active stresses are sufficient ingredients to induce sponta- 
neous biomimetic beating in a clamped active filament. 



FIG. 5. (Color) Logarithm of the absolute value of the char- 
acteristic function F as a function of A and a, obtained ana- 
lytically in the free-draining approximation. The golden lines 
track the zero contours, and correspond to the various normal 
modes as indicated. The transition to instability occurs via a 
saddle-node bifurcation. 
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